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ABSTRACT 


The  effects  of  solar  radiation  heating  on  a  large,  thin,  shallow  spherical  orbiting  shell 
system  are  studied.  The  shell  is  representative  of  a  large  antenna  /reflector  system  which  could 
also  be  gravitationally  stabilized  by  connecting  a  sub-satellite  to  a  long  tether.  After  determining 
the  steady  state  temperature  distribution  across  the  thickness  of  the  shell,  the  corresponding 
thermal  deformation  is  estimated  as  a  function  of  the  solar  incidence  angle  and  material  properties. 
The  closed-loop  dynamic  response  of  the  shell  is  simulated  by  using  LQR  control  design  synthesis 
and  assuming  that  12  point  actuators  are  placed  on  the  outer  surface  and  edge  of  the  shell. 
Calculations  also  indicate  that  the  shell's  transverse  elastic  vibrational  frequencies  may  differ 
significantly  from  their  nominal  values  due  to  the  solar-thermal  influence. 

A  practical  control  strategy  for  the  minimum-time  maneuver  problem  has  been  developed 
and  successfully  applied  to  the  Naval  Research  Laboratory's  Reconfigurable  Spacecraft  Host 
for  Attitude  and  Pointing  Experiments  (RESHAPE)  hardware  test  facility.  This  is  the  first  time 
that  the  bang-bang  type  and  feedback  control  strategy  has  been  applied  to  the  RESHAPE  facility 
by  including  a  flexible  appendage  which  is  represented  with  the  attachment  of  a  spherical 
pendulum  device  to  the  edge  of  the  RESHAPE  platform.  The  tests  resulted  in  an  excellent  • 
correlation  between  the  numerical  simulations  and  the  experimental  test  results,  and  demonstrates 
the  effectiveness  of  the  feed  forward  and  state-error  feedback  control  strategy. 
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INTRODUCTION 


A  tethered  reflector/antenna  system  has  been  suggested  for  possible  use  in  mobile 
communications,  power  transmission  and  astronomical  observations  and  has  been  the  subject  of 
study  in  several  past  AFOSR  grant  reports.  The  tether  would  be  attached  to  a  large  orbiting 
antenna/reflector  and  also  to  a  small  subsatellite.  For  communications  and  power  generation  the 
subsatellite  could  contain  the  electronic  feed  or  power  generating  system.  The  system  could  be 
gravitationally  stabilized  by  deploying  a  small  (200-500  kg.)  subsatellite  a  distance  of  Ikm 
downward  along  the  local  vertical. 

The  second  chapters  of  this  report  extends  the  previous  analyses  of  the  effect  of 
environmental  disturbances  due  to  solar  radiation  pressure  to  now  include  the  heating  attributed 
to  the  thermal  gradients  which  result  in  thermally  induced  deformations.  In  many  cases  there  will 
be  a  noticeable  interaction  between  the  solar  radiation  pressure  and  the  thermally  distorted 
system.  In  this  chapter  the  dynamics  and  control  of  an  orbiting  shallow  spherical  shell  subject  to 
solar-thermal  interaction  is  studied.  Expressions  for  the  thermal  deflections  and  resulting  solar 
radiation  pressure  induced  torques  are  developed.  Control  of  the  orientation  and  shape  of  the 
shell  can  be  accomplished  by  the  use  of  actuators  located  at  certain  points  on  the  shell  surface 
and  outer  edge. 

In  last  year's  grant  report  the  successful  testing  of  three-axes  near  minimum  time 
maneuvers  using  the  Naval  Research  Laboratory's  Reconfigurable  Spacecraft  Host  for  Attitude 
and  Pointing  Experiments  (RESHAPE)  was  reported.  An  excellent  correlation  between 
numerically  calculated  and  experimentally  obtained  results  was  accomplished.  These  previous 
results  were  based  on  a  rigidized  model  of  the  RESHAPE  platform  representative  of  three-axis 
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maneuvers  of  a  rigid  spacecraft  system.  Recently  these  results  haye  been  extended  by  utilizing  a 
pendulum  attached  vertically  by  a  string  to  the  edge  of  the  RESHAPE's  circular  platform  which 
could  represent  a  flexible  appendage.  The  first  degree  of  freedom  of  the  pendulum  can  be 
described  by  a  swing  motion  with  respect  to  a  vertical  line  through  the  attachment  point  in  the 
tangential  plane,  and  represents  the  first  flexible  appendage  mode.  A  second  degree  of  freedom 
is  represented  by  the  out-of-the  tangential  plane  angular  motion.  In  Chapter  Three  of  this  report 
the  recent  numerical  and  experimental  results  of  near-minimum  time  maneuvers  of  the  NRL 
RESHAPE  test  facility  with  the  flexible  pendulum  appendage  are  described.  The  control  strategy 
is  based  on  the  feedforward  (open-loop  solution)  plus  the  state-error  feedback  control;  this 
strategy  has  been  successfully  used  in  the  previous  experimental  study  based  on  the  rigid-model 
of  the  RESHAPE.  This  is  the  first  time  that  a  bang-bang  type  and  feedback  control  strategy  has 
been  applied  to  the  RESHAPE  testing  facility  with  a  flexible  appendage  included. 

Finally  Chapter  Four  summarizes  the  over-all  conclusions  and  suggestions  for  further 

work. 


1.2 


II.  SOLAR  -  THERMAL  EFFECTS  ON  A  SHALLOW  SPHERICAL  SHELL  IN  ORBIT 


Abstract 

The  effects  of  solar  radiation  heating  on  a  large,  shallow,  spherical  shell  in  orbit  are  considered. 
After  determining  the  steady  state  temperature  distribution  across  the  shell  thickness,  the  thermal 
deformation  experienced  is  estimated  as  a  function  of  the  material  properties  and  the  solar 
incidence  angle.  There  is  a  harmonic  variation  of  the  physical  dimensions  of  the  shell  as  a 
function  of  orbital  position.  The  dynamic  response  of  the  shell  is  simulated  by  utilizing  a 
configuration  of  twelve  control  actuators  placed  on  its  outer  surface  and  on  its  edge,  and  utilizing 
the  LQR  design  methodology.  More  significantly,  calculations  show  that  the  frequencies  of 
transverse  vibrations  of  the  shell  elastic  modes  experience  considerable  deviation  fi’om  their 
nominal  values  due  to  solar-thermal  effects.  This  shift  in  the  frequencies  must  be  taken  into 
consideration,  for  large  shells  which  are  inherently  flexible. 

2.1  Introduction 

The  effect  of  environmental  disturbing  influences  on  large  space  structures  in  orbit  has  received  a 
great  deal  of  attention  in  recent  years.  The  major  environmental  disturbances  on  such  structures 
will  be  due  to  solar  radiation  pressure  and  solar  heating.  The  use  of  thin,  shallow,  spherical,  shell 
type  structures  as  principal  components  of  large  space  structures  is  also  quite  widespread.  In 
particular,  shallow  spherical  shell  type  structures  may  be  used  as  reflectors  or  antennas  to 
perform  various  tasks  in  such  disciplines  as  telecommunications,  astronautics,  space  exploration, 
etc.  The  effect  of  solar  radiation  pressure  acting  as  a  load  on  a  class  of  flexible  shallow, 
spherical  shell  type  structures  has  been  considered  previously'.  It  was  determined  that  the  solar 
radiation  pressure  interacts  with  the  elastic  modes  of  the  shell  to  produce  significantly  higher 
rigid  modal  oscillations  as  compared  to  the  oscillations  produced  if  the  shell  were  treated  as 
completely  rigid.  Another  important  environmental  effect  is  the  heat  due  to  solar  radiation  that 
results  in  thermal  gradients  in  a  structure.  For  specific  materials,  the  thermal,  mechanical,  and 
surface  properties  ate  such  that  the  deformations  that  are  caused  by  thermal  gradients  can 
sometimes  result  in  dynamic  instability  of  the  orbiting  structured  In  most  cases,  there  is  also  a 
significant  interaction  between  the  solar  radiation  pressure  and  the  thermally  deformed  structure. 
The  dynamics  of  orbiting  beams  and  plates  in  the  presence  of  such  a  solar-thermal  interaction  has 
also  been  investigated^.  In  this  paper,  the  dynamics  of  an  orbiting  shallow  spherical  shell 
subjected  to  such  an  interaction  will  be  studied.  Expressions  for  the  thermal  deflection  of  a 
spherical  shell  will  be  developed.  The  solar  radiation  pressure  torque  acting  on  such  a  thermally 
deformed  shell  will  be  estimated  and  an  attempt  will  be  made  to  control  the  attitude  of  the  shell 
by  means  of  actuators  strategically  placed  on  its  outer  surface  and  on  the  edge.  Finally,  the 
extent  of  any  thermoelastic  interaction  will  be  investigated  by  considering  the  shift  in  the  natural 
frequencies  of  the  shell  flexible  modes  as  a  result  of  the  solar-thermal  interaction. 
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2.2  Assumptions 

The  general  problem  of  solar-thermal  interaction  and  the  associated  thermoelastic  problem  is 
very  complicated  for  structures  that  span  more  than  one  dimension.  The  problem  becomes  even 
more  complicated  if  thermal  transients,  the  shadowing  of  the  structure  and  the  non-linear 
boundary  conditions  associated  with  radiation  heat-transfer  from  the  exposed  surfaces  are  all 
taken  into  accoimt.  In  the  present  analysis,  the  following  simplifying  assumptions  are  made  in 
order  to  reduce  the  problem  to  manageable  proportions. 


1.  The  thin,  shallow,  spherical  shell  is  assumed  to  move  in  a  100%  sunlit,  circular  orbit  at  a 
geosynchronous  altitude.  Furthermore,  solar  radiation  is  restricted  to  the  plane  of  the 
orbit. 

2.  The  effects  of  the  Earth  reflected  solar  radiation  and  the  thermal  radiation  emanating  from 
the  Earth  are  neglected. 

3 .  The  inherent  time  lags  in  the  heat  transfer  process  are  very  small  compared  with  the 
orbital  period  (=  24  hours  for  a  geosynchronous  orbit)  and  are  ignored. 

4.  Owing  to  the  small  thickness  of  the  shell,  the  radiation  from  its  edge  can  be  neglected. 

5.  The  local  shadowing  of  a  part  of  the  shell  due  to  the  remaining  part  can  be  neglected 
because  of  the  extreme  shallowness  of  the  shell. 

6.  The  steady  state  temperature  distribution  is  strictly  a  function  of  the  radial  coordinate. 
Thus,  the  temperature  varies  only  along  the  shell  thickness.  This  necessarily  implies  that 
at  steady  state,  the  isothermal  surfaces  are  concentric  spheres.  Another  implication  of 
this  assumption  is  that  no  shear  stresses  are  introduced  as  a  result  of  the  thermal  gradient. 
As  a  consequence,  the  thermal  deflection  will  have  a  lone  radial  component. 

7.  For  the  simulation  times  (2-3  orbit  periods)  considered  in  the  analysis,  it  will  be  assumed 
that  the  position  of  the  Sxm  remains  fixed  with  respect  to  an  Earth-fixed  inertial  frame  of 
reference. 


2.3  The  Thermal  Problem 

The  equation  that  governs  the  conduction  of  heat  in  a  spherical  domain  is  given  by: 


dT 

'dt 


kV^  T  =  K  ( • 


dr 


2  dT 


r  dr 


(1) 


where,  k  is  the  thermal  diffusivity,  r  is  the  radial  coordinate  (along  the  shell  thickness),  t  denotes 
time  and  T  =  T  (r,t)  is  the  temperature.  The  steady  state  temperature  distribution,  is  given  by: 


d^T 
r - 

dr^ 


(2) 


The  solution  to  (2)  can  be  achieved'*,  by  changing  the  dependent  variable  to  rT.  If  r^,  r„  and  T,,,  Tj 
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respectively  denote  the  radii  of  curvature  and  the  steady  temperatures  at  the  outer  and  inner 
surfaces  of  the  shell,  then  it  can  be  easily  shown  that: 


T{r) 


r(r^-rj 


(3) 


Fig.  1.  shows  the  cross-sectional  geometry  of  the  shallow  spherical  shell  along  with  it 
dimensions.  Fig.  2.  shows  the  cross-section  of  the  shell  in  the  orbit  plane  in  a  state  of  thermal 
equilibrium.  Here,  t,  denotes  the  unit  vector  in  the  direction  of  the  incoming  solar  flux,  making 
an  instantaneous  angle  of  0;  with  the  outward  normal  to  the  shell.  The  solar  incidence  angle  Oj, 
is  assumed  to  be  a  constant  for  a  small  interval  of  time.  During  this  interval,  the  surface  facing 
the  sun,  S^,  attains  a  steady  state  temperature  T^,  while  the  surface  away  from  the  sxin,  Sj,  attains 
a  steady  state  temperature  Tj.  The  total  heat  leaving  the  shell  from  the  two  surfaces  should  equal 
the  heat  received  by  the  shell  from  the  solar  flux.  At  equilibrium,  the  heat  flowing  through  the 
shell  along  the  radial  direction  must  also  equal  the  heat  radiated  from  the  unexposed  surface  of 
the  shell.  The  above  two  statements  represent  heat  energy  balance  equations,  the  mathematical 
equivalents  of  which  will  be  used  to  solve  for  the  steady  state  temperatures  attained  by  the  outer 
and  the  inner  surfaces  of  the  shell. 


If  %  and  Cj  respectively  denote  the  emissivities  of  the  outer  and  the  inner  surface  of  the  shell,  o, 
the  Stefan-Boltzmann  constant  (=56.7  X  10  '^  KW/m^-K'*),  k,  the  coefficient  of  thermal 
conductivity  of  the  shell  material,  h,  the  nominal  value  for  the  thickness  of  the  shell,  aj,  the 
absorptivity  coefficient  of  the  exposed  shell  surface,  and  G,  the  intensity  of  solar  radiation  (=  0.8 
KW/m^)  for  Earth-orbiting  structures,  then  the  above  mentioned  statements  of  thermal 
equilibrium  result  in: 

4  4 

e  or  +  e^oT.  =  a.  Gcos9,(t) 

o  Q  i  i  5  1 


when  the  outer  surface  of  the  shell  is  exposed  to  solar  radiation  and. 
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When  the  inner  surface  of  the  shell  is  exposed  to  solar  radiation.  In  the  present  analysis, 
depending  on  the  instantaneous  (angular)  location  of  the  shell  in  its  orbit,  the  set  of  equations 
represented  by  either  (4)  or  (5)  will  be  used  to  numerically  solve  for  the  steady  state  temperatures 
attained  by  the  inner  and  outer  surfaces.  The  Newton-Raphson  algorithm  will  be  used  for  this 
purpose.  It  is  clear  from  (4)  and  (5)  that  the  amount  of  thermal  radiation  received  by  the  shell 
from  the  sun  is  greatest  when  the  outward  normal  vector  to  the  exposed  shell  surface  and  the  unit 
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vector  in  the  direction  of  the  incoming  solar  flux  are  parallel  to  each  other.  This  occurs  (when  0; 
=0),  twice  per  orbit.  It  is  reasonable  to  expect  that  the  thermal  gradient,  |  To  -  Tj  |  will  be  greatest 
at  these  two  locations.  The  results  can  be  seen  in  Fig.  3,  which  shows  the  steady  state  surface 
temperatures  as  a  function  of  time,  for  a  shell  made  up  of  polyamide  (k  =  0.245  X  10'^  KW/m- 
®K)  and  when  Cq  =  Cj  =  0.05  and  =  0.2.  It  can  be  seen  that  for  each  orbit,  the  only  period  of 
time  when  the  inner  surface  attains  a  steady  state  temperature  greater  then  the  steady  state 
temperature  of  the  outer  surface  is  when  the  inner  surface  is  exposed  to  the  sun. 

2.4  Thermal  Deformation 


With  the  assumption  that  the  temperature  varies  only  along  the  radial  direction,  it  can  be 
demonstrated^  that  the  thermal  deflection  is  also  a  function  of  the  radial  coordinate  r.  Figure  4, 
shows  an  infinitesimal  element  of  the  shallow  spherical  shell  in  a  state  of  spherically  symmetric 
stress.  It  can  be  demonstrated^  that  under  these  conditions  all  shear  stress  components  vanish 
and  09  =  aj=  which  will  be  referred  to  as  the  tangential  stress  component.  Utilizing  the.  stress 
equilibrium  equations,  the  stress-strain  constitutive  relations  and  the  strain-displacement  field 
equations  in  the  spherical  domain,  it  can  be  shown  thaf : 
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where,  E,  v,  and  a  are  respectively  the  Young's  modulus,  the  Poisson's  ratio  and  the  coefficient 
of  linear  expansion  of  the  material  of  the  shell.  A  and  B  are  constants  to  be  determined  from  the 
boundary  conditions.  Using  (3)  in  (8),  and  making  the  assumption  that  o  does  not  vary  across 
the  shell  thickness,  it  can  be  shown  that: 
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In  order  to  determine  A  and  B  one  can  impose  the  condition  that  the  normal  radial  stress 
components  vanish  at  the  inner  and  outer  surfaces. 

Or(r.)  -  a^(r^)  =  0  (10) 


Using  (3)  in  (6)  results  in, 
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Using  (1 1)  in  (10),  one  can  solve  for  A  and  B  as: 
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Finally,  using  (12)  and  (13)  in  (9)  results  in: 
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Since  the  isothermal  surfaces  inside  the  shell  are  concentric  spheres,  (14),  gives  the  radial 
displacement  of  the  generic  isothermal  at  radius  r  e  [ri,ro].  For  the  problem  at  hand,  the  solar 
radiation  pressure  torques  acting  on  the  thermally  deformed  shell  depend  primarily  on  the  radius 
of  curvature  of  the  inner  and  outer  surfaces.  Thus,  the  radial  displacements  on  these  surfaces  are 
of  special  interest.  It  is  easily  verified  that: 
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The  deformation  given  by  (14)  will  be  along  the  direction  of  the  surface  that  attains  the  greater 
steady  state  temperature.  Because  of  the  symmetric  nature  of  this  deformation,  the  overall  effect 
is  to  cause  an  expansion  of  the  interior  of  the  shell  without  altering  its  shape.  As  the  shell  moves 
in  its  orbit,  the  change  in  the  intensity  of  solar  radiation  to  which  it  is  exposed  causes  it  to 
experience  periodic  heating  and  cooling  phases  which  alternate  once  for  each  orbit  (Fig.3).  As  a 
consequence  of  these  alternate  heating  and  cooling  phases,  the  shell  undergoes  periodic 
expansions  and  contractions  that  cause  the  principal  dimensions  of  the  shell  (a,  q,  r^,  and,  H)  to 
vary,  but  it  retains  its  spherical  shape. 

Fig.  5  shows  the  deformation  undergone  by  the  generic  isothermal  surface  at  r  to  a  position  at  r , 
as  dictated  by  (14).  If  a',  H',  and  i|f'  respectively  denote  the  altered  values  of  the  base  radius,  the 
height,  and  the  negative  of  the  distance  of  the  center  of  the  center  mass  from  its  center  of 
curvature,  then  from  the  geometry  of  Fig.  5,  it  is  obvious  that: 
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Figures  6,7, 8, 9  and  10  respectively  show  the  variation  of  the  shell  inner  radius  of  curvature,  the 
shell  outer  radius  of  curvature,  the  shell  thickness,  the  shell  base  radius  and  the  shell  height  as  a 
function  of  orbital  time.  The  shell  inner  and  outer  radii  experience  significant  variations  with 
time  due  to  solar-thermal  effects.  A  peak  deviation  of  nearly  17  m  from  the  nominal  value  of 
5000  m  is  observed.  Figures  8-10  show  that  the  remaining  dimensions  of  the  shell  only  undergo 
small  changes  in  magnitude. 

2.5  The  Solar  Radiation  Model 

The  solar  radiation  pressure  torques  induced  about  the  center  of  mass  of  the  deformed  shell  will 
form  the  focus  of  attention  of  this  section.  Expressions  for  the  solar  radiation  forces  and  torques 
acting  on  space  structures  of  arbitrary  shape  were  derived  by  Karymov®.  In  the  Karymov 
formulation,  solar  radiation  force  and  torque  expressions  are  first  obtained  for  the  separate  cases 
of  totally  absorbing  and  totally  reflecting  surfaces.  If  N,  and  denote  the  solar  radiation 
torques  on  a  completely  absorbing  and  a  completely  reflecting  surface,  respectively,  then, 

At  =  ii  t  X  fp{t.n)dS 

a  o  J 

s  n  8^ 

At  =  2h  fn  X  p  (t.n)^dS 
r  o  J 

S 


Where,  h^  is  the  solar  radiation  constant  =  4.64  X  10'®  N/m^,  t  is  the  unit  vector  in  the  direction 
of  the  incoming  solar  flux,  p  is  the  position  vector  of  the  generic  point  on  the  shell  surface  from 
its  center  of  mass,  n  is  the  outward  unit  normal  vector  at  the  generic  point,  and  dS  is  the  area  of 
the  infinitesimal  element  on  the  deformed  shell  surface.  The  integration  in  (18)  is  over  the 
surface  of  the  structure  that  is  actually  exposed  to  solar  radiation  and  is  bounded  by  the  condition 
t.n  0.  Fig.  1 1  shows  a  shell  with  outer  radius  of  curvature  R  and  the  coordinate  frame  used  to 
model  solar  radiation  effects.  The  origin  of  the  coordinate  frame  is  at  the  instantaneous  center  of 
curvature  of  the  shell  and  the  z-axis  is  aligned  with  its  symmetry  axis.  From  Fig.  1 1,  the 
following  expressions  can  be  obtained: 

p  =  Ksin<t)cos  01  +  i?  sin  (J)  sin  9/  *  ( Rcos  k 

t  =  -sinacos(3i-sinasin(3  j-cosai? 
n  =  cos  0  sin  01  +  sin0sin0j  +  cos(t>k 
dS  =  K^sin0d0d0 


Where  a  and  p  are  the  two  solar  incidence  angles  that  uniquely  fix  the  orientation  of  t,  and 
(R,0,(j))  are  the  spherical  coordinates  associated  with  a  generic  point  on  the  shell  surface.  Using 
(10)  in  (18),  one  can  derive  the  following  expressions  for  the  solar  radiation  torques: 
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Here,  <|)o  denotes  the  constant  apex  angle  of  the  shell.  Furthermore,  since  the  unit  vector  in  the 
direction  of  the  incoming  solar  flux  is  assumed  to  lie  entirely  within  the  orbit  plane,  p  =  n/2  will 
have  to  be  used  in  (20)  and  (21).  With  this  assumption,  the  solar  radiation  pressure  induces  a 
torque  only  about  the  pitch  axis  of  the  shell.  In  the  Karymov  formulation,  the  overall  solar 
radiation  torque  acting  on  a  shell  whose  surface  has  an  arbitrary  coefficient  of  reflectivity,  e,  is 
given  by: 

W  =  (l-e)  *  eW^  (22) 


As  the  shell  moves  about  in  its  orbit,  the  solar  radiation  torques  induced  (equations  (20)-(21)) 
vary  not  only  because  of  a  variation  in  the  solar  incidence  angles  a  and  p,  but  also  because  of  a 
variation  in  the  dimensions  of  the  shell  due  to  thermal  effects  (Figures  6-10).  Figures  12  and  13 
show  the  variation  of  the  solar  radiation  torque  acting  about  the  pitch  axis  for  the  cases  of  a  shell 
with  a  completely  absorbing  and  a  completely  reflecting  surface,  respectively.  Also  shown  in 
these  figures  is  the  variation  of  the  torques  when  the  thermal  deformation  of  the  shell  is  not 
considered.  The  figures  demonstrate  that  the  inclusion  of  thermal  deformation  effects,  results 
only  in  a  negligible  increase  in  the  solar  radiation  torque  induced. 

2.6  Equations  of  Motion 


The  equations  of  motion  of  an  orbiting  shallow  spherical  shell  were  derived  by  Bainum  and 
Kumar’.  The  linearized  equations  of  a  completely  rigid  shell  can  be  shown  to  be: 
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Here  q^  qy_  and  q^  respectively  denote  the  shell  pitch,  roll,  and  yaw.  J^,  Jy,  and  4  denote  the  mass 
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moments  of  inertia  of  the  shell  about  the  axes  of  the  reference  frame  of  Fig.  1 1 .  C„  Cy,  and 
denote  the  external  torques  acting  about  the  axis  denoted  by  the  subscript.  The  attitude  of  the 
shell  is  assumed  to  be  controlled  by  means  of  12  actuators  placed  on  the  surface  and  on  the  edge 
of  the  shell  in  the  configuration  of  Fig.  14.  In  Fig.  14,  all  actuators  not  on  the  edge  of  the  shell 
produce  forces  directed  along  the  normal  to  the  surface  at  the  point,  whereas  all  actuators  on  the 
edge  produce  forces  that  have  equal  components  along  the  normal  to  the  surface,  and  along  the 
tangent  to  the  edge.  Thus,  the  contribution  to  the  external  torques  C^,  Cy,  in  (23)  is  due  to  the 
solar-thermal  effects  as  well  as  the  restoring  torques  produced  by  the  actuators.  The  equations  of 
motion  given  by  (23)  are  rewritten  in  state-space  form  as  x'  =  AjX  +  B,u,  and  a  full  state  feedback 
control  law  of  the  form  u  =  -Kx,  that  minimizes  the  performance  index: 

J  =  I  X  ^Qx  *  u'^Ru  dt  (24) 

0 


is  sought.  Here,  A,  and  B,  are  the  state  and  the  control  influence  matrices,  respectively,  while  Q 
and  R  are  their  corresponding  weighting  matrices.  Using  Q  =  100  Igxe.  and  R  =  Ii2xi2»  the 
controlled  response  of  the  orbiting  shell  is  simulated  and  the  response  of  the  pitch,  the  roll,  and 
the  yaw  obtained  is  shown  in  Figures  (15),  (16),  and  (17)  respectively. 


2.7  Effect  on  Flexibility 


The  flexibility  of  the  shell  is  modelled  by  assuming  that  the  orbiting  shell  vibrates  primarily 
along  the  transverse  direction.  The  transverse  elastic  vibrations  of  a  thin,  shallow,  spherical  shell 
with  a  completely  free  edge  were  first  derived  by  Reissner  and  Johnson*.  They  derive 
expressions  for  the  natural  frequencies  and  mode  shapes  of  the  shell  and  show  that  each  mode 
shape  is  uniquely  associated  with  a  pair  of  integers  (n,k),  where  n  and  k  refer  to  the  number  of 
nodal  meridians  and  circles,  respectively.  Furthermore,  the  frequency  of  vibration,  On  k,  of  the 
(n,k)th  mode  is  related  to  a  frequency  parameter  ^  according  to  the  relation: 


k 


D 

hpa  ^ 


Ca  ^ 
DR^ 


(25) 


where  h  is  the  shell  thickness,  p  is  the  mass  density,  and  E  is  the  Young's  modulus  of  the  shell 
material.  C  =  Eh  is  the  longitudinal  stiffness  factor,  and  D  =  EhVl2(l  -  v^),  the  bending  stiffness 
factor,  with  v  denoting  the  Poisson's  ratio  of  the  shell  material.  It  can  also  be  demonstrated*,  that 
for  the  first  few  modes  (for  which  the  shell  vibrates  with  zero  or  one  nodal  meridian)  the 
frequency  parameter,  Pn j,  satisfies: 


2 
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which  must  be  numerically  solved  for  k.  Here,  J„(.)  and  !„(.),  respectively  denote  the  Bessel 
function  and  the  modified  Bessel  function  of  the  first  kind  of  order  n.  Equation  (25)  can  then  be 
used  to  solve  for  the  frequency  of  vibration  Wn  k.  In  the  context  of  the  present  problem,  the 
variation  of  the  physical  dimensions  of  the  shell  from  solar-thermal  effects  causes  a  shift  in  the 
natural  frequencies  of  its  elastic  modes.  Figures  (18)-(23)  show  the  variation  of  o)„  ij  with  time 
for  the  first  few  modes.  Also  shown  in  the  figures  is  the  nominal  value  of  for  the  specific 
mode.  For  the  dimensions  of  the  shell  considered,  these  nominal  values  of  the  frequencies  are 
quite  closely  packed.  Consequently,  there  is  a  significant  shift  in  the  frequencies  from  the 
deformation  caused  by  solar-thermal  effects.  In  fact,  for  most  of  the  first  few  modes,  it  is 
observed  that  the  peak  value  of  the  frequency  attained  is  actually  larger  than  the  nominal 
frequency  of  the  next  mode. 

2.8  Conclusions 

In  this  paper,  a  simple  mathematical  model  that  incorporates  the  solar-thermal  effects  on  an 
orbiting  shell  has  been  presented.  It  has  been  demonstrated  that  for  a  shell  made  up  of 
polyamide,  a  peak  thermal  gradient  of  nearly  3.6  °C  is  observed  across  the  1  cm  shell  thickness. 
The  shell  is  alternately  heated  and  cooled  in  its  orbit  in  a  periodic  fashion  due  to  a  periodic 
change  in  the  amount  of  solar  radiation  to  which  it  is  exposed.  The  shell  experiences  a  peak 
increase  in  its  radius  of  curvature  of  nearly  17  m  along  with  smaller  increases  in  its  other 
dimensions.  Although  these  changes  in  the  dimensions  do  not  significantly  increase  the  solar 
radiation  torque  acting  on  the  shell,  they  cause  a  significant  variation  in  the  frequencies  of 
vibration  of  the  shell  flexible  modes.  This  suggests  that  in  the  process  of  the  design  of  control 
laws  of  large  flexible  shells,  particular  emphasis  must  be  placed  in  controlling  the  transverse 
vibrations.  This  is  particularly  true  when  mission  requirements  dictate  that  stringent  pointing 
accuracies  be  satisfied. 
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Shell  Surface  Temperatures  (Degree  Celsius) 


Fig.  1.  Shell  Cross-Sectional  Geometry  Fig.  2.  Thermal  Gradient  Across  Shell  Thickness 

Due  to  Solar  Radiation  Heating 
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Fig.  3.  Variation  of  Steady  State  Fig.  4.  Stress  Distribution  on  an  Infinitesimal 

Surface  Temperatures  Element  in  the  Shell  Interior 


Variation  of  Shell  Inner  Radius  From  Solar-Thermal  Effccts( 36000  Km) 


Fig.  5.  Thermal  Deformation  of  Fig.  6.  Shell  toner  Radius  Vs.  Orbital  Time 

the  generic  Isothermal 
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Shell  Outer  Radius  of  Curvature  (Meters) 


Vanauon  of  Shell  Outer  Radius  From  Solar-Thermal  Errccls(36000  Km) 


0  0.2  0,4  0.6  0.8  I  1.2  1.4  1.6  1.8 

Orbit  Period 

Fig.  7.  Shell  Outer  Radius 
Vs.  Orbital  Time 

Variation  of  Shell  Base  Radius  From  Solar-TIiermal  Effects  (36000  Km) 
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Fig.  9.  Shell  Base  Radius 
Vs.  Orbital  Time 
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Fig.  1 1 .  Reference  Frame  Used  to  Model 


Variation  of  Shell  Tliickness  From  Solar-Thermal  Effects  at 
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Fig.  8.  Shell  Thickness 
Vs,  Orbital  Time 

Virialion  of  Shell  Height  From  Solar-Thermal  Effects  at  36000  Km 
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Fig.  10.  Shell  Height  Vs. 
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Fig.  12.  Solar  Heating  Effects  on  Pitch  Torque 
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Fig.  14.  The  Twelve  Actuator  Configuration 
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m.  LARGE-ANGLE  NEAR-MINIMUM-TIME  MANEUVER  OF 
RESHAPE  WITH  FLEXIBLE  APPENDAGES 


Abstract 

In  this  chapter,  we  present  our  recent  numerical  and  experimental  results  of  the  near¬ 
minimum-time  maneuvers  of  the  Naval  Research  Laboratory's  Reconfigurable  Spacecraft  Host  for 
Attitude  and  Pointing  Experiments  (RESHAPE)  facility  with  flexible  appendages.  These 
maneuvers  are  designed  as  single  axis  maneuvers,  but  are  tested  on  the  three-dimensional  testbed. 
The  control  strategy  is  based  on  the  feedforward  (open-loop  solution)  plus  state-error  feedback 
control.  The  open-loop  solution  is  obtained  by  solving  the  so-called  two-point  boundary-vdue 
problem  associated  with  the  minimum-time  maneuvers  using  a  shooting  method.  This  solution 
provides  a  bang-bang  type  control  command  for  the  prescribed  reorientation  requirements,  i.e.,  to 
slew  the  structure  from  the  known  initial  states  to  the  final  required  attitude  in  a  miiiimum-time 
sense.  In  the  experiment,  this  open-loop  control  solution  provides  a  "feedforward"  control 
command  and  trajectories  which  the  system  is  required  to  follow.  The  second  loop,  a  state-error 
feedback  control  loop,  is  then  designed  to  add  guaranteed  stability  to  the  system,  so  that  small 
disturbances  during  the  maneuver  and  post-maneuver  error  can  be  damped  out.  The  numerical 
and  experimental  results  are  presented  to  illustrate  the  success  of  this  method. 


3.1  Introduction 

The  objective  of  these  maneuver  tests  is  to  achieve  minimum-time  or  near-minimum-time 
maneuvers  of  a  ground-based  space  structure  with  flexible  appendages  and,  at  the  same  time,  to 
keep  the  vibration  of  the  flexible  components  in  a  minimum.  The  Naval  Research  Laboratory's 
RESHAPE  facility  testbed  is  chosen  for  our  test  (  refer  to  Refs.  1-4).  The  RESHAPE's 
configuration  is  shown  in  Fig.  1.  The  testing  platform  is  supported  by  a  spherical  air  bearing 
system  and  the  platform  can  rotate  freely  about  the  pivot  point,  0,  the  center  of  the  sphere.  To 
achieve  the  objective,  a  pendulum  is  attached  vertically  by  a  string  to  the  edge  of  the  RESHAPE's 
circular  platform  to  represent  the  flexible  appendage.  A  sketch  of  the  RESHAPE's  coordinate 
system  with  the  appendage  is  shown  in  Fig.  2.  The  rigid  body  system  is  then  augmented  by  the 
two  degrees  of  freedom  of  the  pendulum  motion  under  the  restoring  force  induced  by  the  gravity 
on  the  pendulum. 

The  first  degree  of  freedom  of  the  pendulum  can  be  characterized  by  the  swing  angle,  a, 
of  the  pendulum  with  respect  to  the  vertical  line  through  the  hinge  point,  a.  The  second  degree  of 
freedom  can  be  represented  by  the  out  of  the  tangential  plane  angular  motion,  i.e.,  by  the  angle 
between  the  projection  line  of  the  pendulum  vector  on  the  platform  plane  and  the  tangential  plane, 
P,  as  shown  in  Fig.  2. 
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In  this  test,  we  consider  only  small  amplitude  flexible  vibration  motion,  since  we  do  not 
want  to  excite  the  flexible  vibration  too  much.  Under  this  consideration,  the  3-dimensional  (3-D) 
equations  of  motion  of  the  system  have  been  derived.  It  can  be  seen  that  the  out-of-plane  equation 
of  motion  is  coupled  with  other  equations  through  a  coefficient  sina,  where  a  is  the  angle  of  the 
pendulum  with  respect  to  the  virtual  line.  Therefore,  in  the  linearization  process  that  follows,  the 
out-of-plane  motion  is  neglected.  As  a  result  we  consider  only  the  swing  motion  of  the  pendulum 
in  the  tangential  plane  which  can  be  effectively  considered  as  one  flexible  mode.  In  our  tests,  we 
mainly  consider  large-angle  maneuvers  about  the  yaw  axis.  The  motion  about  the  roll  axis  and  the 
pitch  axis  are  kept  small.  This  reduces  the  control  design  to  single-axis  maneuver  control  of  a 
rigid  body  with  a  single  flexible  mode.  However,  the  test  is  still  a  3-D  test. 

In  the  open-loop  control  design,  a  shooting  method  is  used  to  solve  the  two-point 
boundary-value  problem  resulted  from  the  Maximum  principle  for  the  maneuver  problem.  This 
solution  provides  a  bang-bang  type  control  command  for  the  prescribed  reorientation 
requirements,  i.e.,  to  slew  the  structure  from  the  known  initial  states  to  the  final  required  attitude 
in  a  minimum-time  sense.  In  the  experiment,  this  open-loop  control  solution  provides  a 
"feedforward"  control  command  and  trajectories  which  the  system  is  required  to  follow.  The 
second  loop,  a  state-error  feedback  control  loop,  is  then  designed  to  add  guaranteed  stability  to 
the  system,  so  that  small  disturbances  during  the  maneuver  and  post-maneuver  error  can  be 
damped  out. 


3.2  System  Equations 

The  3-D  dynamical  equations  of  motion  can  be  derived  using  Lagrange's  approach.  The 
kinetic  energy  of  the  system  can  be  expressed  as 


(1) 


where  the  subscripts  A/,  h,  and  m  refer  to  the  platform,  the  reaction  wheels,  and  the  pendulum, 
respectively.  Then, 


(2) 


where  oo  is  the  angular  velocity  vector  and  is  only  the  inertia  matrix  of  the  platform  about  the 
pivot  point,  O.  Both  (o  and  are  expressed  in  body  axes  (i,  j,  k).  Next,  by  assuming  that  the 
rotation  axes  of  the  three  reaction  wheels  are  aligned  orthogonally  along  the  three  body  axes  (i,  j, 
k),  respectively,  we  can  obtain 
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where  is  the  inertia  matrix  of  the  /th  wheel  with  respect  to  the  pivot  point ,  O,  Jj  is  the  moment 
of  inertia  of  the  /th  wheel  about  its  rotation  axis;  and  ^  is  the  relative  angular  rate  of  the  /th 
wheel.  The  kinetic  energy  related  to  the  pendulum  can  be  obtained  as 


^X+2(oX^+Z)i+— (oV^to 
m 


where  m  is  the  mass  of  the  pendulum,  J„  is  the  inertia  matrix  of  the  pendulum  with  respect  to  O, 
and  ( referring  Fig.  2) 
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where  Rj  is  the  magnitude  of  the  vector  R  projected  onto  the  (i,  j)  plane  of  the  RESHAPE 
coordinate  system  (Fig.  2).  The  "d"  sign  represents  a  3x3  skew-symmetric  matrix  associated  with 
a  vector  a. 
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The  virtual  work  related  to  the  gravity  can  be  obtained  as 


hW=-W[(M*m^gh  +  mgiR^L)\a^-mg{bttLl+h^Ll)a^  (5) 


where  60  is  the  virtual  angle  displacement  vector,  is  the  third  column  of  the  direction  cosine 
matrix  of  the  platform. 


A  = 


2(?i4'2-?o^3)  2(?o+<?2Vi  2(^2^3+4rp^i) 
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M  is  the  mass  of  the  platform,  is  the  mass  of  three  wheels,  b  is  the  position  vector  from  the 
pivot  point  to  the  center  of  mass  of  the  system  without  the  pendulum,  g  is  the  gravitational 
acceleration  constant,  and  q  is  the  attitude  quaternion  vector,  q^={  qo  qj  q:  ^3  ]■ 

We  now  can  obtain  Lagrange's  equations  for  the  quasi-coordinates  related  to  co,  the 
rotational  angles,  (p,  of  the  reaction  wheels,  and  the  pendulum  swing  angles,  a  and  P, 


+  m(R+L)L  = 


dco 
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where 
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is  the  torque  acting  on  the  reaction  wheels,  Jm+hIS  the  inertia  matrix  of  the  platform  and  the 
wheels  together,  and  Qg  is  the  generalized  force  vector.  By  substituting  the  generalized  forces 
from  Eq.  (5)  and  Eq.  (7)  into  Eq.  (6),  we  can  remove  the  term  Q, 


=  U+&  — — m(L)-j^a>-  [(M +tnj)gb  +mg(R  *L)\a^ 
d(i> 


(8) 


At  the  same  time,  we  obtain  the  equations  for  a  and  P, 


i  1 

m[Z^Z  -Zj(A+Z)(b-2ZjZ<o-— &)— co]  =  -mga^L^ 
“  m  da 


(9) 


T*  1  T 

w[Z  ZQ-Z»(/i+Z)<b-2Z»Z<o — <0— «]  =  -mga^L^ 

m  dp 


(10) 


cosasinP 

sinacosP 

z  = 

a 

cosacosP 

-smasmP 

sina 

0 

It  can  be  proven  that 
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It  can  be  seen  here  that  for  small  angle  a  and  P, 


cosa-1,  cosP=cos(  p-p^)“COsP|^ 
sina-a,  smp=sm(p-p^)“-smp^ 


the  P  motion  Eq.  (10)  is  coupled  Avith  w  motion  Eq.  (8)  and  a  motion  Eq.  (9)  by  a  small  quantity 
sina.  Therefore,  for  our  present  small  a  amplitude  analysis,  we  drop  Eq.  (10)  from  further 
consideration. 

For  the  purpose  of  simulation,  we  need  to  obtain  the  linearized  equation  of  the  system.  In 
addition  to  small  a  and  P,  we  also  assume  small  q, ,  q2,  Wj  and  coj.  As  a  result,  all  second  and 
higher  order  terms  related  to  a,  P,  qi,  q2,  ,  and  (O2  are  neglected. 
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Note  that  these  are  linear  equations  for  a,  P,  qi,  q2,  Wi ,  and  (O2  only.  Jy  is  the  element  of  the 
constant  inertia  matrix  Jj- , 

where 
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and 


l=L^-R^,  c=XjCosP^,  ir=I,smP^ 

To  further  simplify  the  analysis,  we  can  neglect,  the  terms  related  to  the  gravity,  the  body  rates, 
and  the  rates  of  the  wheels  in  the  o)  equation,  and  the  terms  related  to  the  body  rates  and 
quaternions  in  the  a  equation. 
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The  2-dimensional  (2-D)  equations  involving  only  (O3  and  a  can  be  obtained  immediately  as. 
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mJJjZjCbj  +  miiO  =  -wgZjO 


(15) 


(16) 


3.3  Time  Optimal  Control  of  the  2-D  System  and  3-D  Feedback  Control 

The  time  optimal  control  for  the  2-D  system  can  be  obtained  by  using  the  maximum 
principle.  The  formulation  of  the  time-optimal  control  is  based  on  the  following  cost  functional 
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The  Hamiltonian  of  the  system  can  be  formed  by  introducing  costates  and  the  necessary 
conditions  for  the  optimal  control  can  be  derived.  The  control  is  of  bang-bang  type  and  can  be 
obtained  numerically  by  solving  the  associated  two-point  boundary-value  problem.  The  shooting 
method  from  Ref  2  is  used  for  solving  the  nonsingular  TPBVP. 

The  solution  of  the  optimal  control  of  the  2-D  system  is  used  as  the  open-loop  control  or 
feedforward  control  of  the  3-D  system.  This  solution  includes  the  quaternion,  q(t),  angular 
velocity,  co(t),  and  control,  u(t).  The  solution  for  a(t)  is  not  considered  in  the  feedback  control 
loop.  The  errors  between  these  computed  values  and  the  measured  values  ("real")  of  these 
variables  are 


«<o  =  -  <o(0 


where  (»)„  is  obtained  from  the  output  of  the  rate  gyros,  and  q„,  is  obtained  from  integration  of 
the  kinematic  equations  using  Euler's  first  order  integration  method.  These  errors,  the  error 
quaternion,  6q,  and  error  angular  velocity,  6  a),  are  used  in  the  state-error  feedback  control  loop 
to  obtain  the  correction  control,  6u, 

f  =  l,  2,  3. 


where  kp  and  kp  are  constants.  The  resulting  control  signal  is  the  combination  of  the  open-loop 
control  and  the  feedback  control. 


«(r,  bq,  6 to)  =  u(t)  +  bu 


It  is  clear  that  the  torque  constraints  for  the  open  loop  control  (bang-bang  type)  should  be  set 
below  the  maximum  allowable  level  of  the  physical  system,  so  that  some  leeway  can  be  left  for  the 
correction  control.  As  a  rule  of  thumb,  we  can  use  80-90%  of  the  allowable  control  for  the  open- 
loop  control  design  and  10-20%  for  the  correction  control.  This  results  in  the  near-minimum-time 
maneuvers  under  the  allowable  control  constraints. 
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3.4  System  Parameters 


The  description  of  the  RESHAPE  facility  can  be  found  from  Ref  1.  The  RESHAPE 
structural  parameters  used  in  the  simulations  are  listed  below: 

Inertial  matrix  of  the  RESHAPE  without  the  pendulum: 


14.42 

.15 

.28 

.15 

17.85 

.42 

.28 

.42 

27.82 

(ft-lb-Sti) 


Momentum  of  inertia  of  the  reaction  wheels:  Ji=  J2  =  J3  -.054  (ft-lb-sec^ ); 

Gravitational  constant:  g  =  32.174  (ft/sec^); 

RESHAPE  weight  =  638.0  (lb); 

Angle  between  the  position  vector  of  the  pendulum  hinge  point  and  body  i  axis:  P=7r/2; 

Length  of  the  position  vector  of  the  pendulum  hinge  point:  Ri=2.5  (ft); 
z-coordinate  of  the  position  vector  of  the  pendulum:  R^  =  0.0; 

Length  of  the  pendulum:  Li  =  4.5  (ft); 

Pendulum  weight:  mg  =  5.0  (lb); 

Position  vector  of  the  RESHAPE  without  the  pendulum:  bi=0.0;  bj  =  -mRi/  (  M  +  m^ );  b3=0.0; 
Sampling  Time:  0.05  sec; 

Feedback  control  gains:  kp  =[24.8846  30.2138  49.7636],  k,  =[44.22  53.69  88.43]; 
Maximum  control  torque  for  the  open-loop  control  Uj :  =  0.4  (ft-lb); 

Extra  control  effort  for  the  feedback  control  are: 

6U3  =  0.001056-0.08976  (ft-lb),  depending  on  the  cases; 

Extra  control  for  the  other  two  axes  (pitch  and  roll)  are:  6  U2  =  6  Uj  =  0.56  U3 


3.5  2-D  Optimal  Control  Solution  and  3-D  Simulation 

The  2-D  dynamic  equations  (15)  and  (16)  are  used  to  find  the  time-optimal  control 
solution.  This  is  a  single-axis  maneuver  of  a  rigid  body  with  a  flexible  mode.  We  here  consider 
only  the  rest-to-rest  maneuvers  and  the  associated  boundary  conditions  of  the  state  variables  are 
chosen  as 


(0,(0)  =  0,  63(0)  =  630; 

6>3(9  =  0,  63(9  =  0. 


a(0)  =  0,  d(0)  =  0, 

a(9  =  0,  &(9  =  0, 


3.9 


For  the  RESHAPE  system  with  parameters  described  in  section  3.4,  we  have  obtained  maneuver 
solutions  from  830  =  1  deg  to  830  =  180  deg.  Two  typical  maneuver  solutions,  25  deg  and  180  deg 
maneuvers,  are  shown  in  Fig.  3  and  Fig.  4,  respectively.  The  time-optimal  controls  for  these 
maneuvers  are  of  bang-bang  type  and  three  intermediate  switches  are  required.  The  minimum  time 
for  the  25  deg  maneuver  is  1 1.312  sec  and  three  intermediate  switches  are  separated  by  0.6  sec. 
The  minimum  time  for  the  180  deg  maneuver  is  30.091  sec  and  the  time  interval  between  the 
three  intermediate  switches  is  0.6  sec.  The  occurrence  of  triple  switches  is  clearly  associated  with 
the  addition  of  the  flexible  mode  included  into  the  system  compared  with  the  purely  single-axis 
maneuver  of  a  rigid  body.  The  well-known  result  for  the  latter  case  is  a  single  switch  in  the  middle 
of  the  slewing  period.  It  is  shown  that  the  pendulum  angle  and  pendulum  rate  are  very  small  for 
this  particular  RESHAPE  system  setup.  This  small  amplitude  vibration  is  also  observed  during  the 
hardware  test. 

We  have  observed,  through  computations,  that  the  time  interval  between  the  triple 
switches  is  different  from  case  to  case,  depending  on  the  system  parameters  and  initial  conditions. 
If  all  other  parameters  are  flxed  and  only  the  initial  condition,  830,  is  changed  from  1  deg  to  180 
deg,  the  time  interval  between  the  three  switches  changes  periodically.  At  some  particular  values 
of  830,  we  are  able  to  obtain  a  zero  interval.  This  means  a  single  switch  case  where  all  the  three 
switches  are  collapsed  together,  even  though  there  is  a  flexible  mode  involved  here.  But  these  are 
only  a  few  isolated  cases  and  an  elaborate  analysis  will  not  be  conducted  in  this  report.  Except 
these  isolated  points,  all  other  cases  result  in  triple  switches.  We  choose  the  25  deg  maneuver 
because  the  triple  switches  in  this  case  are  well  separated,  and  choose  the  180  deg  case  because  it 
can  be  considered  as  the  "largest"  angle  rotation  for  the  rest-to-rest  maneuvers. 

The  solution  for  these  maneuvers  is  usually  not  available  in  analytical  form  compared  with 
the  classical  single-axis  rigid  body  maneuver,  and  therefore,  can  not  be  generated  at  run-time. 
However,  we  are  interested  in  finding  out  the  advantage  of  applying  this  solution  to  the  maneuver 
test  over  using  the  single  switch  solution.  Before  we  go  to  the  test,  let  us  look  at  results  of  the  3- 
D  simulations  for  these  two  open-loop  control  commands.  The  results  for  the  triple  switch  and  the 
single  switch  results  are  presented  in  Figs.  5,  and  6  for  the  25  deg  maneuver  case,  respectively. 

The  single  switch  command  for  the  latter  method  is  obtained  as  if  there  is  no  flexible  part.  By 
comparing  these  two  figures,  we  can  see  that  using  the  triple  switch  command  will  reduce  the 
vibration  of  the  pendulum  both  during  the  second  half  of  the  maneuver  and  post  maneuver.  This  in 
return  will  reduce  the  attitude  errors  and  rates.  The  pendulum  angle  and  rate  are  also  shown  in 
Figs.  5  and  6.  The  testing  results  for  this  comparison  will  be  discussed  later. 

We  should  point  out  that  we  use  both  control  command  feedforward  as  well  as  state 
trajectory  feedforward  in  our  control  strategy.  This  is  quite  different  from  using  only  the 
command  feedforward  method.  Actually,  we  can  observe  this  by  examining  Fig.  5  and  Fig.  6. 
Before  the  first  switch  near  5  sec.,  both  control  commands  (open-loop)  are  the  same  U3=-0.4  fl-lb. 
But  the  3-D  responses  of  the  control,  U3  are  different.  Also  the  6q3  in  these  two  figures  are 
different.  This  is  because  the  feedforward  trajectories  for  these  two  method  are  different,  as 
shown  in  Fig.  7. 
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The  time  responses  of  the  reaction  wheels  for  the  triple  switch  simulation,  (p  and  Q  are 
also  shown  in  Fig.  7  for  the  case  6u=0.03696.  Other  cases  are  similar.  The  nonzero  values  for  Q 
are  due  to  the  gravitational  torque  terms  presented  in  Eq.  (1 1).  This  phenomenon  is  also  observed 
during  the  test,  where  the  reaction  wheels  for  the  x  and  y  axes  keep  running  after  the  maneuver.  If 
the  values  of  a  or  q  are  not  zero,  there  will  be  some  gravitational  torques  on  the  platform  which 
balance  the  torques  from  the  reaction  wheels.  Even  though  the  torques  applied  to  the  wheels  are 
zero,  other  parameters  or  factors  may  also  lead  to  the  nonzero  momentum  accumulation  of  the 
wheels  during  the  maneuver.  Further  analysis  is  needed  on  this  subject. 


3.6  3-D  Test  Results  and  Numerical  Simulations 

Using  the  feedforward  and  state-error  feedback  method  described  in  section  3.3,  we  have 
conducted  the  hardware  tests  and  3-D  numerical  simulations.  Five  tests  for  the  25  deg  maneuver 
are  presented  here  as  shown  in  Figs.  8-11.  The  first  four  tests  are  using  the  triple  switch  solution 
and  the  fifth  test  is  using  the  single  switch  method  for  the  comparison  purpose.  For  the  triple 
switch  maneuvers,  we  have  used  four  different  extra  control  effort  levels,  6U3,  for  the  feedback 
control  use  as  discussed  in  section  3.3.  These  values  are  0.001056,  0.01056,  0.03696,  0.08976 
(ft-lb)  for  the  four  cases,  about  0.2%  to  22%  of  the  open-loop  control  saturation  level,  0.4  (ft-lb). 
Also  shown  in  Figs.  8-1 1  are  the  3-D  simulations  using  the  3-D  equations  (7),  (1 1),  and  (12). 

From  Figs.  8-1 1,  as  the  extra  control  effort  is  increased,  the  errors  of  are  reduced  and 
the  overshoot  of  (O3  at  the  end  of  the  maneuver  is  also  reduced.  The  errors  for  qj,  qj,  &>i,  and  CO2 
do  not  change  significantly.  We  also  see  that  there  is  no  similarity  between  testing  results  and 
simulation  results  for  errors  of  qi,  qj,  Wj,  W3,  Uj,  and  U2.  Further  analysis  needs  to  be  done  in  this 
area.  We  also  see  that  the  control  history  of  the  simulation  is  close  to  the  open-loop  control 
command  (not  shown),  while  the  control  signal  from  the  test  deviates  more  from  the  open-loop 
solution  either  for  the  smaller  leeway  case.  Fig.  8c,  or  the  larger  leeway  case.  Fig.  1  Ic.  We  believe 
that  the  overshoot  of  the  control  at  the  end  of  the  maneuver  is  due  to  the  inaccuracy  of  the  inertial 
term  and  the  difference  in  the  inertial  terms  between  the  2-D  and  3-D  representing  models  of  the 
system. 

In  Fig.  12,  we  put  the  triple  switch  case  and  single  switch  results  together  for  comparison. 
We  observe  that  the  triple  switch  method  results  in  less  residual  vibration,  even  when  the  error 
magnitudes  are  the  same.  This  indicates  the  advantage  of  using  the  triple  switch  solution 
especially  for  the  well-separated  time  interval  between  the  three  middle  switches. 

The  test  for  180  deg  maneuver  are  presented  here  in  Figs.  13-15.  Three  different  levels  of 
6U3  are  considered.  Again,  either  small  leeway  or  large  leeway  results  in  large  deviation  of  the 
control  from  the  bang-bang  type.  The  case,  6u3=0.03696  (ft-lb,  about  9.2%  of  0.4  ft-lb),  in  Fig. 

14  seems  to  be  the  closest  one  to  the  bang-bang  type  control  of  the  open-loop  solution.  It  is  also 
seen  that  the  errors  in  the  quaternion  during  the  maneuver  are  reduced  as  more  control  freedom  is 
allowed. 
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Another  problem  with  the  180  deg  maneuver  is  that  the  triple  switches  presented  in  the 
simulations  do  not  occur  in  the  test.  The  reason  for  this  needs  to  be  found  out.  One  reason  might 
be  that  the  180  deg  maneuver  requires  longer  slewing  time  than  the  25  deg  maneuver.  The  "mid" 
point  or  middle  region  is  "relatively"  small  compared  with  the  whole  slewing  time  and  is, 
therefore,  not  significant  for  the  feedback  control  adjustment. 


3.7  Conclusions 

A  practical  control  strategy  for  the  minimum-time  maneuver  problem  has  been  developed 
and  successfully  applied  to  the  RESHAPE  hardware  test  facility.  The  open-loop  time-optimal 
control  provides  a  basis  for  the  nonlinear  large-angle  near-minimUm-time  maneuver,  while  the 
linear  feedback  control  takes  care  of  disturbances  during  the  maneuver  and  damps  out  the  residual 
vibration  after  the  maneuver.  This  is  the  first  time  that  the  bang-bang  type  and  feedback  control 
strategy  has  been  applied  to  the  RESHAPE  hardware  testing  facility  with  flexible  appendages. 

An  extensive  optimal  control  solution  for  the  2-D  system  has  been  obtained.  3-D 
simulations  by  using  both  triple  and  single  switch  methods  have  been  conducted.  Use  of  the  triple 
switch  command  will  reduce  the  vibration  of  the  pendulum  and  the  attitude  errors  both  during  the 
second  half  of  the  maneuver  and  post  maneuver.  As  the  extra  control  effort  is  increased,  the  state 
errors  are  reduced.  An  adequate  level  of  the  extra  control  leeway  seems  to  be  about  1 0%  of  the 
open-loop  control  saturation  level  so  that  the  control  will  be  close  to  the  bang-bang  type,  the 
maneuver  will  be  near  minimum-time,  and  less  post  maneuver  vibration  and  state  errors  will 
result. 


The  success  of  the  test  is  demonstrated  by  the  excellent  correlation  between  numerical 
simulations  and  experimental  testing  results.  This  test  demonstrates  the  applicability  and  the 
effectiveness  of  this  simple  feedforward  and  state-error  feedback  control  strategy.  This  technique 
can  be  easily  transferred  to  other  types  of  three-axis,  nonlinear,  and  time-optimal  maneuver 
experiments  or  actual  tests. 

Further  analysis  needs  to  be  done  in  the  areas  such  as  the  identification  of  the  isolated 
single  switching  point,  the  reduction  of  the  terminal  angular  rates  of  the  reaction  wheels,  the 
determination  of  the  reason  for  the  non-similarity  of  the  state  errors  between  the  test  and  the 
simulations,  and  the  reason  for  the  disappearance  of  the  triple  switches  in  the  test  for  the  180  deg 
maneuver  case.  Improvements  and  extensions  of  the  control  strategy  to  three-axis  slewings  and 
vibration  suppression  of  the  rigid  platform  containing  other  flexible  appendages  or  antennas  are 
recommended. 
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Notebook  PC 


Fig.  1  Configuration  of  the  RESHAPE  Facility,  NRL 
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(a)  Vector  Notations.  (b)  Top  View. 
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(c)  Side  View. 

Fig.  2  RESHAPE  Coordinate  System. 
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Fig.  7  3-D  Maneuver  Simulation;  25  deg,  du3=0.03696  ft-lb; 

Feedforward  Trajectory  Difference  between  the  Triple  and  Single  Switch  Cases; 
Reaction  Wheel  Responses  for  the  Triple  Switch  Case. 
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Fig.  9a 
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Fig.  10a 

RESHAPE  Maneuver: 

25  deg  Simulation  and  Test 
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Fig.  10b 

RESHAPE  Maneuver: 

25  deg  Simulation  and  Test 

du=0.03696  ft-lb 
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Fig.  12  RESHAPE  Maneuver  25  deg  Tests: 

Single  Switch  Case:  du3>0. 08976  ft-lb; 
Triple  Switch  Case:  du3=0.08976  ft-lb; 
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Fig.  13  RESHAPE  Maneuver  180 


Test,  Time  (s) 


Simulation  and  Test:  du3=0.01056  ft-lb. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  effects  of  solar  induced  heating  on  a  large,  thin,  shallow  spherical  shell  in  orbit  are 
modelled  mathematically.  Such  a  structure  would  be  representative  of  proposed  orbiting 
antenna/reflector  systems.  For  a  100  m.  base  radius  shell  made  of  polyamide  a  maximum  thermal 
gradient  of  3.6“ C  is  predicted  across  the  1  cm  shell  thickness.  As  the  system  moves  in  its  orbit 
the  shell  is  alternately  heated  and  cooled  in  a  periodic  rhanner  due  to  the  periodic  change  in  the 
amount  of  solar  radiation  to  which  it  is  exposed.  The  shell  is  predicted  to  experience  a 
maximum  increase  in  its  radius  of  curvature  of  17m  (nominal  value,  5000  m)  with  smaller 
increases  in  its  other  dimensions.  There  is  no  corresponding  significant  increase  in  the  solar 
radiation  torque  acting  on  the  shell  but  there  is  a  significant  change  in  the  shell's  flexible 
vibrational  frequencies.  In  fact  for  most  of  the  first  few  modes  the  peak  value  of  the  frequency 
attained  in  the  presence  of  solar-thermal  effects  is  actually  larger  than  the  nominal  frequency  of 
the  next  highest  mode.  This  means  that  when  the  mission  requirements  stipulate  very  accurate 
pointing  and  shape  control  that  particular  emphasis  should  be  placed  in  the  design  of  the 
controllers  that  limit  the  transverse  vibrations. 

Possible  extensions  to  this  work  could  consider  the  thermal  shock  problem  encountered 
when  the  system  transfers  from  the  sunlit  side  of  the  orbit  to  the  shadowed  side  and  also  the  effect 
of  local  shadowing  by  different  subsystem  components  on  the  rest  of  the  system  elements. 
Another  extension  to  this  work  is  to  develop  a  model  that  could  predict  the  solar-thermal  effects 
on  the  tethered-shallow  spherical  shell-subsatellite  system.  This  would  represent  an  extension  to 
the  research  reported  in  Chapter  III  of  Volume  I  of  this  report,  which  emphasized  the  solar 
radiation  effects  on  the  tethered  reflector  sub-satellite  system. 
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A  practical  bang-bang  type  and  feedback  control  strategy  for  the  minimum-time  maneuver 
problem  has  been  successfully  applied  for  the  first  time  to  the  Naval  Research  Laboratory's 
RESHAPE  hardware  test  facility  with  a  flexible  pendulum-type  appendage.  The  control  strategy 
is  based  on  the  feed  forward  (open-loop  solution)  plus  state-error  feedback  control.  The  open- 
loop  solution  is  obtained  by  solving  the  related  two-point  boundary  value  problem  using  a 
shooting  method.  The  success  of  this  testing  is  shown  by  the  excellent  correlation  between  the 
numerical  simulations  and  the  experimental  test  results. 

The  technique  developed  here  can  be  easily  applied  to  other  types  of  three-axis,  nonlinear, 
and  time-optimal  maneuver  experiements.  It  is  suggested  that  this  technique  could  be  applied  to 
more  sophisticated  hardware  facilities  such  as  the  Phillips  Laboratory's  ASTREX  faciility  in  the 
future.  Unfortunately,  during  this  AFOSR  grant  year  it  was  not  possible  to  schedule  testing  time 
with  the  original  ASTREX  facility  at  Edwards  AFB  although  preliminary  data  was  exchanged 
and  numerical  simulations  were  performed.  Perhaps,  with  the  relocation  of  the  ASTREX  faciility 
to  Albuquerque,  NM,  another  opportunity  might  be  available.  In  addition  to  possible  future 
testing  with  the  ASTREX,  three-axis  slewing  and  vibration  suppressionof  the  RESHAPE  platform 
with  other  types  of  flexible  appendages  such  as  beams  or  antennas  is  recommended.  Depending 
on  the  type  of  attachment  of  the  appendages  the  effect  of  gravity  may  have  to  be  compensated 
such  as  with  the  use  of  follower  support  cables  and  wires.  Such  devices  would  then  have  to  be 
included  in  the  RESHAPE  system  mathematical  models,  similar  to  the  way  that  they  are  currently 
included  in  the  ASTREX  models. 
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